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In  the  density  transition  layer  between  salt  and  fresh  water  high 
frequency  interfacial  waves  can  be  generated  by  the  shear  stresses  created 
by  finite  amplitude  internal  waves  at  this  interface.  The  conditions 
under  which  these  waves  occur  are  experimentally  investigated  in  this 
thesis.  Instantaneous  Richardson  number  profiles  at  the  crests  and 
troughs  of  the  finite  amplitude  internal  waves  were  measured.  The  low 
frequency  waves  were  mechanically  generated  by  a hydraulically  driven 
horizontal  diaphragm. 

The  experiments  were  performed  in  an  internal  wave  facility  with  a 
test  section  eighty  feet  long,  four  feet  high  and  two  feet  wide. 
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CHAPTER  I 
INTRODUCTION 

Mass  density  strati f ication  in  a fluid  plays  an  essential  role  in  the 
gravity  wave  phenomenon.  As  stated  by  C.  S.  Yih  (1965),  "wave  motion  can 
only  exist  in  a stably  stratified  fluid".  Indeed,  a free  surface  is  just 
an  extreme  case  of  density  stratification  with  a density  discontinuity 
at  the  surface. 

From  a physical  point  of  view,  the  oscillatory  wave  motion  is  caused 
by  the  restoring  force  acting  on  a fluid  particle  as  the  particle  is  dis- 
placed from  its  moan  position.  The  restoring  force  which  causes  the 
particle  to  seek  its  way  to  its  mean  position  is  the  buoyancy  force.  So 
buoyancy  forces  are  the  primary  forces  associated  with  internal  gravity 
waves  that  appear  at  the  interface  between  two  fluid  layers. 

Internal  waves  are  a common  and  important  phenomenon  in  the  ocean  as 
well  as  in  the  atmosphere.  They  can  have  a significant  effect  on  both 
fixed  and  floating  submarine  structures.  Internal  waves  are  thought  to 
be  a major  source  for  the  production  of  turbulence  in  the  oceans  and 
atmosphere.  Theoretical  investigations  have  shown  that  internal  waves 
can  be  generated  by  a number  of  difference  mechanisms.  Under  certain  con- 
ditions shear  flows  in  a stratified  fluid  can  be  unstable,  and  the  growth 
of  this  instability  can  result  in  internal  waves.  Once  finite  amplitude 
internal  waves  exist,  they  in  turn  can  generate  higher  frequency  internal 
waves  as  follows. 

In  internal  wave  motion  the  fluid  particle  velocity  above  the  inter- 
face is  in  the  opposite  directii  n to  the  particle  velocity  below  the  in- 
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terface,  thus  a large  velocity  gradient  occurs  at  the  interface.  This 
can  provide  the  necessary  conditions  for  a disturbance  to  grow,  i.e.  an 
unstable  region  in  the  flow.  J.  D.  Woods  (1973)  observed  such  wave  in- 
duced disturbances  in  the  ocean  and  referred  to  them  as  wave-induced  shear 
instabilities.  This  is  similar  to  those  observed  by  Thorpe  ( 1 909 ) in 
the  laboratory,  and  the  clear-air  turbulence  in  the  atmosphere  identified 
by  Gossard,  et.  al . (1970). 

In  previous  experimental  studies  of  the  mixing  due  to  internal  waves, 
by  0.  M.  Sheppard,  T.  Toniasello,  and  l.  Chou,  high  frequency  disturbances 
(or  interfacial  waves)  were  observed  at  the  interface.  The  disturbed 
region  was  confined  to  a thin  layer  with  clearly  defined  boundaries. 

This  shear  flow  is  unsteady  and  the  interfacial  waves  that  have  been 
observed  are  finite  in  amplitude.  There  has  been  very  little  theoretical 
work  done  on  non-linear  stability  in  general , and  to  the  knowledge  of  the 
author,  no  theoretical  investigation  done  on  non-linear  stability  of  unsteady 
stratified  shear  flows  The  objective  of  this  thesis  is  then  to  obtain 
information  about  the  physical  processes  that  take  place  in  this  complex 
flow  using  the  results  of  the  linear  stability  theory  as  a guide. 
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CHAPTER  II 
BACKGROUND 


Hydrodynamic  stability  has  become  a large  branch  of  fluid  dynamics. 
Stability  of  parallel  flow  of  inviscid  fluid  was  first  studied  by  Helm- 
holtz, Kelvin  and  Rayleigh. 

There  are  several  different  modes  of  instability.  Brooks  Benjamin 
(1963)  has  classified  unstable  differences  in  flexible  surfaces  bounding 
inviscid  flows  into  three  different  types  as  follows: 

(1).  Class  "A"  instability  results  in  Tol lmien-Schl ichting  waves . 

The  presence  of  tin's  type  instability  causes  a dissipation  of 
energy  of  (In'  system  in  am  amount  proportional  to  the  amplitude 
of  the  waves  squared. 

(?).  Class  "R"  instability  is  identified  with  free  surface  waves. 

The  presence  of  this  instability  causes  a raise  in  energy 
level,  as  in  the  case  of  wind  generated  waves. 

(3).  Class  "C"  is  the  well  known  Kel vin-Helmhol tz  instability.  It. 
is  a consequence  of  conservative  forces  acting  on  a small  dis- 
turbance. It  leads  to  a violent  breakdown  at  the  interface. 

The  stability  problem  is  most  readily  treated  in  terms  of  the  theory 
of  small  amplitude  perturbations.  The  linearized  stability  equation  for 
parallel  flows  of  homogeneous  fluids  was  first  derived  by  Orr  and  Sommer- 
fold. This  equation  can  be  written  as 

(u-cHr-k-*)  - I."*  - (*""  - 2k?*  " -t  k"i  ), 
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where  u is  the  primary  flow  velocity,  c the  complex  wave  speed,  <$>  the 
cross-flow  dependence  of  the  disturbance,  k the  wave  number  and  R the 
Reynolds  number.  The  primes  denote  derivatives  with  respect  to  the  cross- 
flow  coordinate. 

The  linearized  governing  equation  for  the  stability  of  inviscid,  den- 
sity stratified  shear  flows  was  first  derived  by  Taylor  and  Goldstein  and 
can  be  wri tten  as 

- [(u-c)  u"  + (u-c)  (u‘)  R.j  + k ] <j>  = 0. 

Here  the  stability  parameter  is  the  gradient  Richardson  number,  R^ , 
which  is  defined  as 

r = . _JLS_  . 

1 p(u')2 

The  mass  density  is  denoted  by  p. 

A sufficient  condition  for  stability  according  to  Miles  and  Howard 
(1966),  is  that  R^  > 1/4.  It  should  be  emphasized  that  this  criterion 
applies  to  infinitesimal  disturbances  in  a parallel  stratified  shear  flow 
of  an  inviscid  fluid.  Woods,  in  a series  of  field  investigations,  found 
that  transition  from  turbulent  to  laminar  flow  occurred  at  Richardson  num- 
bers of  approximately  unity. 

A number  of  analytical  studies  have  been  made  using  hypothetical  ve- 
locity and  density  profiles.  Examples  include  Howard  and  Drazin’s  sin- 
usoidal flow  in  which  the  velocity  profile  is  a sine  function;  Taylor  and 
Goldstein's  shear  layer  in  which  there  is  a linear  variation  of  velocity 
with  the  vertical  coordinate;  and  Holmboe's  discontinuous  density  profile 
in  which  the  density  profile  is  a step  function  at  the  interface.  Hazel 
(1972)  did  some  numerical  studies  of  the  stability  of  stratified  shear 
flows.  He  used  hyperbolic  tangent  functions  for  the  density  and  velocity 
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profiles.  The  error  function  shear  layer  has  also  been  studied.  Most 
of  these  results  are  presented  as  neutral  curves  in  the  k-J  plane,  where 
k is  the  wave  number  of  the  disturbance  and  J is  the  overall  Richardson 
number  based  on  the  total  density  and  velocity  changes 


j . 9(»Wp)h 

(au)2 


A neutral  curve  is  a line  which  separates  the  stable  region  from  the 
unstable  region.  Some  of  the  results  of  these  theoretical  investi- 
gations are  given  in  the  appendix  (figures  64  and  65). 

The  gradient  Richardson  number  can  be  written  as 


*,  ■ - N'V  . 

(■’■“)■ 
a Z ' 

where  N is  the  so  called  Brunt-Vaisala  or  buoyancy  frequency  defined  by 


N = < - ,l/?  • 


This  frequency  plays  an  important  role  in  the  stability  of  stratified 
fluids.  Note  that  the  gradient  Richardson  number  is  a ratio  of  buoyancy 
forces  to  inertia  forces. 

Phillips  (1969)  derived  the  equation 

3‘  ( ) + N2  ( z ) 3‘-w  = 0 

3 t?  3 X2  3 Z2  3 Z2 

from  the  momentum  equations. 

for  small  disturbances  the  vertical  velocity  w can  be  found  by  taking 
w - W( z)  exp  { i ( k x - n t ) 1 

where  k is  the  horizontal  wave  number.  Then  it  follows 

d W + ( N (z)  j-  n-  | kpw  0 
dz2  n2 
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If  n>N,  W is  a monotonic  function  of  z,  the  only  modes  are  surface  waves. 
If  ii'N,  then 


W w exp  (im/) 


where 

. , N-  - n?  «i/? 
m = k ( ) 

o 

n* 

is  the  vertical  wave  number.  W is  an  oscillatory  function  of  z,  it 
corresponds  to  internal  waves,  from  the  above  equation, 


n = N ( — ^ )'/?  . 

k?  + m? 

The  denominator  (k'  » m')’/-1  is  the  resultant  of  the  wave  number  com- 
ponents in  the  x and  z direction  (see  Figure  1). 

(m?  + k7)'/? 


k 

Figure  1 

So  the  frequency  of  the  internal  wave  n = Ncoso  where  |o|  v 1/2  it . 

This  shows  that  when  n >-  N,  internal  waves  travel  at  an  angle  0 from 
the  horizontal  direction,  until  they  reach  a level  where  n N. 

Thus  the  buoyancy  frequency  N is  the  upper  limit  of  the  frequency  for 
which  internal  wave  motions  can  exist  in  a stratified  fluid.  In  a fluid 
with  variable  N,  waves  will  be  trapped  within  the  layer  where  N exceeds 
the  imposed  frequency. 

The  passage  of  an  internal  wave  along  a sharp  interface  between  stra- 
tified layers  produces  a vortex  sheet  at  the  interface.  When  the  wave 
period  is  long  compared  to  the  time  needed  for  any  disturbance  to  grow. 
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an  instability  of  the  Kelvin-Helmhol tz  (K-H)  type  is  to  be  expected. 

If  the  energy  supply  to  the  internal  wave  mode  continues  after  the 
critical  minimum  R.  is  attained,  a local  instability  may  develop.  The 
energy  acquired  by  these  disturbances  is  at  the  expense  of  the  internal 
wave,  so  that  the  occurrence  of  an  instability  of  this  kind  can  restrict 
the  amplitude  of  the  internal  wave. 

In  general  the  breakdown  has  two  effects,  one  destabilizing  and  one 
stabilizing.  If  the  breakdown  is  severe,  mixing  can  occur  which  reduces  the 
density  gradient  thus  destabilizes  the  flow.  The  velocity  gradient  is  like- 
wise reduced  but  this  has  a stabilizing  effect.  It  should  be  pointed  out 
that  stability  parameter  is  more  sensitive  to  changes  in  velocity  gradient 
than  density  gradient. 

Woods  (1969)  observed  the  same  laminar-turbulent-laminar  transition 
in  the  ocean.  He  also  represented  schematically  the  evolution  of  the 
spontaneous  K-H  billows  in  a thermocline  sheet  and  their  subsequent 
transition  to  turbulence.  The  density  and  velocity  differences  across 
the  sheet  did  not  change  during  the  life  time  of  the  disturbances.  The 
mechanism  could  be  described  as  follows. 

After  the  K-H  billow  formed  and  broke,  turbulence  in  the  sheet 
entrained  water  from  above  and  below  thickening  the  sheet.  But  the 
turbulent  sheet  died  when  the  Richardson  number  of  the  sheet  approached 
unity.  Then  the  current  shear  above  and  below  the  sheet  entrains  water 
from  the  thickened  but  no  longer  turbulent  sheet  so  that  it  slowly  thins 
and  reverts  to  laminar  flow. 

Yih-Ho  Pao  (1969)  measured  the  spectra  of  internal  waves  and  tur- 
bulence in  a stably  stratified  atmosphere  and  ocean.  He  classified  in- 
ternal waves  into  3 ranges  according  to  the  wave  number. 
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1).  At  low  wave  numbers,  buoyancy  forces  dominate  over  inertial 
forces  and  the  disturbance  is  in  the  form  of  internal  waves. 

Pao  called  this  the  internal  subrange. 

?).  At  intermediate  wave  number,  the  buoyancy  effect  is  strong  but. 
inertial  forces  are  also  important  and  turbulence  can  exist.  He 
called  this  the  buoyancy  subrange. 

3).  At  high  wave  numbers  the  inertial  forces  are  dominant.  He  called 
this  the  inertial  and  viscous  subrange. 

Although  a great  number  of  analytic  studies  have  been  made,  few  exper- 
imental investigations  have  been  conducted.  Scotti  and  Corcos  (19/1) 
created  a stably  stratified  shear  layer  within  the  test  section  of  a wind 
tunnel  by  merging  two  streams  of  air  after  uniformly  heating  the  top 
stream.  The  two  streams  were  accelerated  side  by  side.  Small  periodic 
disturbances  were  introduced  upstream  of  the  test  section  by  a fine  wire 
oscillating  in  the  thermocl ine.  The  critical  R^  was  found  as  0.2?.  The 
growth  rate  is  a function  of  wave  number.  Wang  (Id??)  first  used  fresh 
water  and  dyed  salt  water  to  form  two  layers  of  fluid  in  a tank.  Shear 
flow  was  induced  by  pumping  the  fluid  in  both  layers.  The  velocity  in  the 
lower  layer  was  larger.  Small  disturbances  were  generated  by  a vibrating 
ribbon  at  the  interface.  Wang  compared  his  results  with  Hazel's  numerical 
analysis  which  showed  good  confirmation.  Prowand  and  Winant  (197?)  con- 
ducted the  same  experiment  in  a two  layer  system,  and  they  observed  the 
growth  of  the  disturbance.  At  t ho  upstream,  the  Richardson  number  is  small 
(R^  = 0.0?  at  x = 0.394  in)  and  the  waves  collapsed.  The  Richardson  rv  - 
her  is  approximately  0.4  at  a distance  of  23.6  in.  Further  downstream, 
the  waves  decayed  and  the  shear  layer  approached  an  undisturbed  laminar 


flow. 


CHAPTER  III 

DESCRIPTION  OF  THE  PROBLEM 


The  objective  of  this  thesis  is  to  investigate  the  conditions  under 
which  high  frequency  interfacial  waves  will  occur  in  the  interfacial  layer 
of  longer  internal  waves.  Long  waves  were  generated  mechanically  in  a 
two  layer  system.  The  upper  layer  was  fresh  water  and  the  lower  layer  was 
salt  water,  dyed  with  red  food  color. 

The  interface  is  defined  as  the  level  where  the  density  gradient  is 
maximum  in  the  two-layer  system.  The  wave  induced  velocity  above  and 
below  the  interface  is  in  opposite  directions  thus  causing  a high  ve- 
locity gradient  (and  shear  stress)  at  the  interface.  According  to  the 
linear,  inviscid  stability  theory,  the  Richardson  number  (R^),  as  de- 
fined earlier,  is  the  stability  parameter  for  an  inviscid,  stratified 
shear  flow.  From  the  definition  of  R^ , a large  density  stratification 
has  a stabilizing  influence  on  the  flow  and  a large  velocity  gradient 
has  a destabilizing  influence. 

The  velocity  gradient  is  maximum  at  the  wave  crests  and  troughs  so 
these  are  points  of  minimum  stability.  In  order  to  determine  if  the 
high  frequency  waves  can  be  generated  by  a shear  instability,  Richardson 
number  profiles  at  the  crests  and  troughs  are  needed.  Velocity,  conduc- 
tivity, and  temperature  profiles  must  then  be  measured  at  these  points. 

Energy  is  transferred  from  the  long  waves  to  the  interfacial  waves 
and  this  results  in  a reduction  in  amplitude  of  the  long  waves.  The  long 
wave  amplitude  is  monitored  during  the  experiments  to  see  if  this  transfer 
of  energy  can  be  correlated  with  the  stability  parameter. 
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A flow  definition  sketch  is  given  in  Figure  ?. 

Two  coordinate  systems  are  used,  one  (/.  the  vertical  coordinate)  has 
its  origin  fixed  to  the  bottom  of  the  tank  and  the  other  (•.,  the  vertical 
coordinate)  has  its  origin  attached  to  the  interface.  Both  i and  s are 
positive  upward. 


CHAPTER  IV 

FACILITIES  AND  INSTRUMENTATION 
Eacjl ities 

Test  Section 

The  experiments  were  conducted  in  the  internal  wave  tank  in  the 
Coastal  and  Oceanographic  Engineering  Laboratory,  at  the  University  of 
Florida.  The  test  section  of  the  tank  is  80  feet  long,  6 feet  high  and 
2 feet  wide  in  the  1 ower  portion  and  3 feet  wide  in  the  upper  wind  tunnel 
section  (see  Figures  3 and  4).  The  tank  is  lined  on  both  sides  with 
1/2  inch  glass  panels  to  facilitate  visual  observation.  The  entire  wave 
tank  is  elevated  18  inches  above  the  ground  to  provide  bottom  access  to 
the  test  section.  The  tank  is  capable  of  generating  wind  waves  and  shear 
flows  and  mechanically  generating  surface  and  internal  wave s. 

Wa ve_Ge n er a tor 

There  is  an  internal  v/ave  generator  at  the  north  end  of  the  wave 
tank  (see  Figure  5).  It  is  comprised  of  a horizontal  38  inch  by  38  inch 
neoprene  rubber  disphragm  attached  to  the  walls  at  its  outer  perimeter, 
and  is  activated  by  a vertical  shaft  at  the  center.  The  shaft  is  driven 
hydraul ically.  A hydraulic  power  supply  is  used  in  conjunction  with  the 
hydraulic  actuator  and  associated  electronic  equipment  to  drive  the  dia- 
phragm. The  frequency  and  amplitude  of  the  waves  are  adjusted  by  a Low 
Frequency  Signal  Generator  Model  202A,  manufactured  by  Hewlett  Packard. 
Storage  Ta nks 

Two  4000  gallon  storage  tanks  are  located  outside  the  building.  These 
tanks  are  connected  to  the  wave  tank  by  4 inch  PVC  pipe.  The  storage 
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Figure  4.  Cross-section  of  the  Internal  Wave  Tank 
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tanks  are  elevated  so  tho  wave  tank  can  60  fill  oil  (com  those  tanks  by 
gravity.  A gear  pump  powered  by  a 6 horsepower  electric  motor  is  used 
to  pump  the  water  back  into  the  storage  tanks. 

Tow  lank  tor  Hot -I  i 1 in  I’rohe  Calibration 

A 11?  inch  by  1.1  inch  by  1?  inch  tow  tank  is  used  to  calibrate  the 
hot-tilm  probe.  The  tank  is  made  of  plexiglass.  A 1/8  horsepower  var- 
iable speed  motor  is  used  to  drive  the  carriage  which  holds  the  hot-film 
probe . 

Traverse  Mechanism 

Itie  traversing  median  i mii  is  comprised  ol  a 1/?  horsepower  variable 
speed  O.C.  motor,  a gear  reduction  box,  a stainless  steel  tube  with  a 
gear  rack  and  a motor  controller  system.  Ihe  traverse  speed  can  be 
varied  from  aprroximately  ;ero  to  6 in/sec.  The  direction  of  the  tra- 
verse shaft  can  be  switched  manually  or  can  be  made  to  shuttle  between 
any  two  positions.  Ihe  velocity  probe,  temperature  probe,  and  conduc- 
tivity probe  were  mounted  at  the  same  level  on  the  bottom  of  the  traverse 
tube. 


Inst rumentat ion 

Position  Indicator 

Ihe  traverse  drives  a ten  turn,  10  k potentiometer  when  it  moves. 
The  potentiometer  is  part  of  a O.C.  circuit  designed  to  give  a linear 
output  voltage  proportional  to  the  position  ot  the  traverse  (see  figure 
66  in  the  appendix).  Ihe  vertical  position  can  be  measured  accurately 
to  0.0?f>  inch. 

Voloc i tv  Measurements 

Velocity  of  the  fluid  is  measured  by  a Constant  Temperature  Ho t - 
I i 1 m Anemometer,  manufactured  bv  1 lw  rmo -Systems , Inc..  Ihe  anomometrs 
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system  consisted  of  a Model  1051  Monitor  and  Power  Supply  module,  two 
Model  1050  Anemometer  modules,  and  two  Model  1055  Linearizers.  The  probe 
which  was  mounted  on  the  tip  of  the  traverse  was  a wedge-type  Model  1233 
NaCl  hot-tilm  sensor.  I he  hot  film  probe  was  calibrated  in  a tow  tank. 
Calibration  curves  were  attained  for  several  salt  concentrations  (specific 
gravities  ranging  from  1.00  to  1.03).  The  change  in  the  probe  output  over 
this  range  of  concentrations  was  found  to  be  loss  than  21'..  Problems 
with  bubble  formation  on  the  probe  at  lower  velocities  in  the  salt  sol- 
utions wore  overcome  by  lowering  the  over-heat  ratio  to  1.05.  The  cal- 
ibration curves  are  shown  in  the  appendix  (Figure  07  and  68  ) . 

1 emperatiirr  Measurement  s 

Temperature  was  measured  with  a glass  encapsulated  thermistor,  ex- 
perimental Type  08127401,  manufactured  by  Keystone  Carbon  Company.  A 
constant  temperature  bath  with  a thermometer  that  could  be  read  to  0.1°F 
was  used  to  calibrate  the  thermistor.  The  calibration  curve  is  shown  in 
the  appendix  (Figure  60) . 

Pens i I \ Measurement s 

Density  strati f icat ion  in  these  experiments  was  achieved  by  varying 
the  salt  concentration.  Since  the  electrical  conductivity  and  mass  den- 
sity of  a salt  solution  varies  with  salt  concentration,  the  mass  density 
of  the  solution  can  be  determined  from  conductivity  measurements.  The 
conductivity  measuring  system  used  in  these  experiments  was  developed 
for  this  work  by  Prs.  Sheppard  and  Poddington  (submitted  for  publication 
in  the  Review  of  Scientific  Instruments).  It  consists  of  a two  elec- 
trode probe  (see  ( igure  b and  7)  and  a constant  peak-to-peak  current  sup- 
ply. Ihe  fluid  being  measured  is  drawn  into  the  probe  through  the  hori- 
zontal slit  at  the  bottom,  through  the  0.012  inch  diameter  hole  and  into 


Figure  6.  Cross-section  of  the  Conductivity  I’rohc 
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the  inner  cavity.  The  fluid  then  goes  out  the  top  of  the  probe  through 
plastic  tubing  to  the  pump.  One  electrode  is  located  in  the  inner  cav- 
ity and  the  other  on  the  exterior  of  the  probe.  The  electrical  current 
density  between  the  electrodes  is  small  except  in  the  0.01?  inch  diameter 
vertical  hole,  where  it  is  quite  large.  Thus  the  conductivity  of  the 
fluid  in  this  hole  at  any  instant  determines  the  resistance  between  the 
el  ec.t.rodos . 

Calibration  of  the  Conductivity  Probe 

To  minimize  error,  calibration  is  done  in  the  tank  under  the  conditions 
present  during  a test.  A siphoning  probe  with  a horizontal  slot  very  sim- 
ilar to  that  in  the  conductivity  probe  is  mounted  on  the  traverse  shaft 
at  the  same  elevation  as  the  other  probes  (conductivity,  thermistor  and 
hot-film).  The  tank  is  then  filled  in  such  a manner  that  t ho  interface 
between  the  salt  and  fresh  water  is  very  diffused.  At  several  points  in 
the  interfacial  layer,  fluid  is  siphoned  from  the  tank  and  conductivity 
readings  recorded.  A hydrometer  (with  0.0001  divisions)  is  used  to  mea- 
sure Hu'  specific  gravity  and  a mercury  thermometer  to  measure  the  temp- 
erature. The  temperature  in  the  tank  is  measured  with  the  thermistor.  A 
calibration  curve  for  the  range  of  salinities  of  interest  can  then  be 
plotted  (see  figure  70  in  the  appendix). 

Internal  Wave  Gages 

The  wave  gages  used  in  this  experiment  were  designed  to  monitor  the 
motion  of  the  interface  between  the  fresh  and  dyed  salt  water  layers. 

The  measuring  system  consists  of  a line  voltage  regulator,  a collimated 
light  source,  a custom  built  light  sensor,  a Wheatstone  bridge  with  a O.C. 
power  supply  and  a recorder.  The  back  edge  of  the  light  sensor  is  para- 
bolic in  shape  (see  figure  8 and  o).  All  surfaces  on  the  sensor 


g ure  3.  Drawing  of  the  Internal  Wave  Gage 
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except  the  front  and  hack  are  painted  black.  The  parabolic  back  surface 
is  polished  and  coated  with  a highly  reflective  paint  and  the  front  light 
admitting  surface  is  just  polished.  A cadmium  sulfide  photo  resistor  is 
placed  on  the  focal  point  of  the  parabola.  The  light  source  which  is  a 
vertical  flourescent  lamp  is  centered  at  the  undisturbed  interface  and 
projects  through  the  wave  tank  to  the  sensor  as  shown  in  Figure  9. 

Honeycomb  is  placed  in  front  of  the  lamps  to  make  the  light  parallel. 

The  salt  water  is  dyed  with  red  food  color  to  take  advantage  of 
the  low  sensitivity  of  the  photo  sensor  to  red  light.  As  the  interface 
moves  up  and  down,  the  amount  of  light  striking  the  light  sensor  changes 
proportionately. 

In  these  experiments,  the  wave  gages  are  only  used  to  measure  the 
wave  frequency  and  wave  length. 

Visicordor 

All  signals  were  recorded  on  a Visicordor  Oscillograph  Model  90GA, 
manufactured  by  Honeywell  as  well  as  on  a magnetic  tape  recorder.  It  has 
a built-in  Time  Line  System  which  provides  accurately  spaced  pulses  at 
intervals  of  0.01,  0.1  and  1.0  second.  It  contains  17  channels  and 
records  on  photographic  paper  with  a high  intensity  ultraviolet  light 
reflected  by  the  galvanometer  mirror. 

Tape  Recorder 

A four-channel,  FM  tape  recorder.  Model  3960A,  manufactured  by  Hewlett 
Packard  was  used  to  record  the  output  voltage  from  the  position  indicator, 
thermistor,  hot-film  probe  and  conductivity  probe. 

The  Pre-Amplifier  and  Low-Pass  Filter 

A four-channel  pre-amplifier  and  low-pass  filter  is  placed  in  series  with 
the  four  sensors.  This  amplifier  is  designed  to  filter  out  electrical  noise 
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induced  hy  equipment  in  the  laboratory.  The  amplifiers  were  used  to 
increase  the  probe  output  to  the  maximum  allowable  input  voltage  of  the 
tape  recorder.  This  unit  has  a high  input-impedance  and  low  output- 
impedance  which  isolates  the  sensors  from  the  recording  equipment.  The 
filter  is  fourth-order,  Butterworth,  low-pass  filter  with  a cutoff  fre- 
quency of  16.7  Hz.  The  circuit  is  shown  schematically  in  the  Appendix 
(Figure  71).  The  low-pass  filter  was  not  used  in  the  hot-film  probe 
circuit. 

X-Y  Plotter 

An  Omnigraphic  ?000  Recorder,  type  6,  manufactured  by  Houston  In- 
struments was  used  during  the  tests  to  plot  the  conductivity  vs.  position. 
This  was  used  to  locate  the  probe  relative  to  the  interface. 


i 


CHAPTFR  V 
PROCEDURE 


Initially  the  wave  tank  was  filled  with  24  inches  of  fresh  water  from 
the  fresh  water  storage  tank.  Salt  was  dissolved  in  the  salt  water  stor- 
age tank  until  the  specific  gravity  of  the  fluid  was  about  1.03.  Red  food 
color  was  used  to  dye  the  salt  water,  and  compressed  air  was  used  to  mix 
the  salt,  dye  and  water.  After  the  fresh  water  was  filled,  the  dyed  salt 
water  was  allowed  to  flow  into  the  tank  very  slowly  from  the  bottom.  After 
6 to  7 hours,  the  salt  water  reached  the  desired  depth  of  24  inches  and 
the  valve  was  closed.  This  slow  filling  procedure  resulted  in  a sharp  in- 
terface between  the  fluids.  After  the  conductivity  probe  was  calibrated, 
the  experimentation  was  begun.  The  conductivity  probe,  thermistor,  and 
hot-film  probe  were  mounted  at  the  same  elevation  on  the  traverse  shaft 
thus  allowing  continuous  simultaneous  measurement  of  salinity,  temperature, 
velocity  and  position.  The  measurement  system  is  shown  in  figure  10.  The 
signals  from  the  probe  were  fed  into  the  4-channel  ampl if ier-f i 1 ter.  The 
filter  can  be  switched  out  of  the  circuit  and  was  not  used  for  the  velocity 
signal.  The  output  signals  were  fed  to  a 4-channel  FM  tape  recorder  and 
along  with  two  wave  gage  outputs  to  a six-channel  pre-amplifier.  The  out- 
puts from  the  pre-amplifier  were  fed  into  an  optical  strip  chart  recorder 
(visicorder) . An  X-Y  plotter  was  used  to  plot  conductivity  vs.  position 
during  the  test. 

An  initial  traverse  was  made  to  obtain  the  initial  density  and  temp- 
erature profiles.  '« he  mechanical  wave  generator  was  then  started  and 
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waves  of  the  desired  frequency  and  amplitude  were  generated.  Since  in- 
stantaneous profiles  could  not  be  made  with  our  measuring  system,  "average" 
profiles  were  measured  using  the  following  procedures.  The  probes  were  po- 
sitioned at  one  point  while  approximately  3 to  5 waves  passed,  then  lowered 
to  a new  position  for  the  next  3 to  5 waves.  This  procedure  v/as  repeated 
until  the  vertical  field  of  interest  had  been  covered.  The  signal  out- 
puts varied  periodically  with  the  wave  motion.  The  point  of  maximum 
S.G.  indicated  the  S.G.  at  the  wave  crest,  and  the  minimum  S.G.  indicated 
the  S.G.  at  the  wave  trough.  So  the  velocity  at  the  crest  was  obtained 
by  reading  the  velocity  at  the  instant  when  the  specific  gravity  was  the 
maximum  and  similarily  the  trough  velocities  were  taken  at  the  point 
where  the  S.G.  was  a minimum.  Since  the  wedge  type  hot-film  probe  used 
in  these  experiments  can  only  measure  velocities  in  one  direction,  only 
velocities  in  the  direction  of  the  propagating  waves  could  be  accurately 
obtained.  Thus  only  profiles  below  the  crests  and  above  the  troughs 
could  be  measured. 

Specific  gravity  profiles  were  constructed  from  the  averaged  values 
at  each  point;  i.e.  the  average  of  the  values  measured  for  3 to  5 con- 
secutive waves. 

Due  to  a high  noise  level  on  the  magnetic  tape  (attributed  to  the 
use  of  old  tape)  the  various  outputs  could  not  be  plotted  on  the  X-Y 
plotter  as  was  originally  planned.  The  visicorder  (optical  strip  chart 
recorder)  output  had  to  be  used,  making  the  data  reduction  a tedious  pro- 
cess. 

The  output  from  two  of  the  optical  internal  wave  gages  were  recorded. 
Knowing  the  spacing  between  these  gages  and  the  approximate  wave  length 
of  the  waves,  the  celerity  and  exact  wave  lengths  were  obtained. 
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Several  amplitudes  wore  examined  for  each  frequency.  This  was  accomp- 
lished without  stopping  the  wave  generator.  The  generator  was  stopped, 
however,  when  the  frequency  was  changed  and  a continuous  traverse  made  at 
that  time. 

Attempts  were  made  to  record  t ho  frequencies  of  the  trapped  waves  at 
the  interface  hut  this  information  was  lost  due  to  the  noisy  tape.  A 
spectral  analysis  to  determine  the  amplitudes  of  the  frequencies  present 
had  been  planned  hut  now  must  he  postponed. 

A table  listing  the  sequence  of  the  various  runs  and  the  duration 
of  each  test  is  given  in  Table  1.  Throughout  the  experiment,  visual 
observations  were  made  and  recorded  and  photographs  were  taken  of 
shadow-graphs  of  the  interface. 
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TABLE  1 


CHAPTER  VI 

RESULTS  AND  DISCUSSION 


Two  typos  of  results  are  presented  here,  visual  observation  and  ex- 
perimental measurements.  A description  of  the  visual  observations  will 
be  presented  first. 

It  was  observed  that  mcchnical ly-generatcd  internal  waves  with 
small  steepness  ratios  (H/L)  maintained  their  smooth  sinusoidal  shape 
both  in  time  and  space.  However,  as  the  steepness  ratio  was  increased 
beyond  some  value,  high  frequency  waves  appeared  in  the  interfacial  layer. 
Upon  the  formation  of  these  interfacial  waves  the  amplitude  of  the  internal 
wave  was  reduced.  After  a period  of  time  the  interfacial  waves  ceased 
to  exist  and  the  internal  wave  amplitude  increased  once  again.  A short 
time  later  the  interfacial  waves  reappeared  and  the  cycle  was  repeated, 
figure  11  shows  sketches  of  this  observed  phenomenon. 

These  high  frequency  waves  were  confined  to  the  interfacial  layer,  i.e. 
tae  transition  layer  for  the  density,  by  very  distinct  boundaries  (see  fig- 
ure 14).  The  interface  was  completely  filled  with  waves  which  from  time  to 
time  would  roll  up,  break  and  appear  to  cause  substantial  mixing  in  the  in- 
terfacial  layer  (see  I inures  13  and  14).  The  measurements  proved  this  not 
to  be  true,  however,  since  only  minor  changes  in  the  density  profiles  were 
observed  even  after  several  hours  of  performing  the  experiments. 

To  aid  in  the  discussion  of  this  phenomenon,  the  following  terminology 
will  be  used.  If  no  interfacial  waves  occur  the  long  internal  waves  are 
said  to  be  stable.  It  interfacial  waves  occur  but  only  after  a long  period 
of  time  this  will  be  referred  t • as  the  transition  conditions.  Waves 


No  interracial  waves.  Large  amplitude 


Interfacial  waves  occur.  Reduced  amplitude 


More  interfacial  waves.  Small  amplitude 


Interfacial  waves  reduced.  Increased  amplitude 


Interfacial  waves  disappear,  Large  amplitude  Velocity  Profile 
figure  11.  Sketch  of  Periodical  Occurrence  of  Interfacial  Waves 
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that  have  high  frequency  internal  waves  appearing  at  the  interface  after 
a short  time  will  be  called  unstable  waves.  The  observed  results  are 
summarized  in  Table  2. 

TA1U.T  2 


Run  No.  Wave  Length 

Wave  Height 

Wave  Period  Wave  Steepness 

Remark 

L (in.) 

llt  (in.) 

T (sec) 

H/L  (10-?) 

2 

117.7 

1 .05 

10.15 

0.8936 

Stabl e 

4 

111.2 

1 .8 

10.15 

1.6187 

6 

108.0 

2.75 

10.15 

2.5463 

Unstable 

9 

52.47 

0.7 

6.85 

1.334 

Stable 

11 

52.5 

1 .2 

6.85 

2.785 

Transition 

15 

40.06 

0.35 

5.23 

0.8737 

17 

28.62 

0.85 

5.23 

2.0369 

Unstable 

The  wave  heights 

, H( , presented 

in  Table 

2 were  measured  after 

allowing 

sufficient  time  for 

the  in  ter facial 

waves  to 

f orm. 

The  experimental 

data  is  presented  next. 

Complete  conductivity  and 

temperature  prof i los 

were  measured 

with  the  wave  maker  turned  off 

before 

and  after  each  sequence  of  experiments.  The  results  of  these  traverses 
which  are  denoted  as  run  numbers  1,  5,  8,  14  and  19  are  presented  in 
Figures  15  through  21.  lor  comparison  purposes  these  profiles  are  plotted 
together  in  Figures  19  and  22. 

The  velocity  anil  density  profiles  were  measured  at  the  wave  crests  and 
troughs.  The  results  are  presented  in  Figures  23  through  43.  Velocity 
profiles  at  wave  crests  were  obtained  below  the  interface  only,  and  pro- 
files at  wave  troughs  were  obtained  above  the  interface  only.  Some  of  the 
velocity  profiles  were  lost  due  to  the  chart  paper  jamming  in  the  visi- 
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Figure  23.  Density  and  Velocity  Profiles  at.  Wave  Crest,  Run  it 2 
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T ho  wave  height  H,  can  be  determined  by  superimposing  the  density 
profiles  at  the  wave  crests  and  troughs.  It  could  also  be  obtained  from 
the  conductivity  position  plots  made  by  the  X-Y  plotter  during  the  ex- 
periments (see  Figure  72  in  the  Appendix).  The  wave  heights  H.  are  com- 
pared with  the  initial  wave  heights  H in  Table  3. 


TABLE  3 


Run  No. 

FI 

0 

(in) 

"t 

(in) 

AM 

(in) 

Remark 

2 

1 .06 

1 .05 

0.01 

Stabl  c 

4 

2.25 

1.8 

0.45 

Transition 

6 

3.25 

2.75 

0.50 

Unstable 

g 

0.87 

0.7 

0.17 

Stable 

11 

1.38 

1.2 

0.18 

Transition 

15 

0.63 

0.35 

0.28 

Trans i tion 

17 

1 .25 

0.85 

0.40 

Unstabl e 

A simple  computer  program  was  written  to  calculate  the  Richardson 
number.  Richardson  number  and  Prunt-Vaisala  frequency  profiles  are  pre- 
sented in  Figures  44  through  Figure  68.  The  minimum  value  on  the 
Richardson  number  profile  for  each  experiment  is  tabulated  in  Table  4. 
These  values  will  be  denoted  by  R^.  . It  should  be  noted  that  in  Figures 
44  through  Figure  58,  the  Richardson  number  becomes  zero  at  points  above 
and  below  the  interface.  This  results  from  the  density  gradient  going 
to  zero  at  these  points.  For  a homogeneous  fluid  the  Richardson  number 
has  no  meaning  and  thus  ceases  to  be  a stability  parameter. 

Interfacial  waves  were  not  generated  by  waves  with  small  steepness 
ratios.  As  the  steepness  ratio  was  increased  the  interfacial  waves  formed 
and  the  internal  wave  height  was  reduced.  The  energy  in  internal  waves 
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is  proportional  to  the  square  of  the  wave  height.  Since  the  wave  maker 
constantly  supplies  energy  to  these  waves,  a reduction  in  wave  height 
would  indicate  that  energy  is  being  dissipated  or  transferred  to  another 
part  of  the  flow,  from  Table  3 it  ran  be  seen  that  every  reduction  in 
wave  height  was  preceded  by  the  formation  of  interfacial  waves.  It  thus 
appears  that  energy  is  transferred  from  the  long  internal  waves  to  the 
interfacial  waves  once  unstable  conditions  are  reached. 

TABl f 4 


•d 

Run  No. 

Min  R. 
i 

Remark 

o 

C 

15.34 

Stahl e 

4 

1.03 

Trans i tion 

6 

1 .96 

On stahl e 

11 

5.55 

Transi t ion 

15 

3.19 

Trans i tion 

17 

2 . 34 

Unstahl e 

When  the  wave  height 

was  reduced 

the  steepness 

velocity  gradient  in  the  interfacial  layer  wore  reduced  accordingly.  This 


i 


I 


would  have  a stabilizing  effect  on  the  flow.  As  would  be  expected,  the 
interfacial  waves  ceased  to  exist  after  a period  of  time  (approximately 
5 to  10  minutes)  and  the  internal  wave  height  increased.  The  stabilizing 
effect  of  the  reduction  in  steepness  ratio  caused  t tie  interfacial  waves 
to  transfer  their  energy  back  to  the  internal  waves.  There  was  very  lit- 
tle energy  dissipated  in  turbulence  as  evidenced  by  the  essentially  con- 
stant density  profile  in  the  interfacial  layer  throughout  the  experiment 
( see  F igure  19) . 
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With  the  increased  steepness  the  flow  becomes  unstable  once  again 

and  the  cycle  is  repeated.  The  energy  transferred  to  the  interfacial 

? ? 

waves  is  proportional  to  H - H . Tin's  energy  transfer  is  plotted 
versus  minimum  Richardson  number  in  Figure  59.  It  appears  that 

? O 

II  - li^  increases  exponentially  with  the  decreasing  minimum  Richardson 
number.  Thus  the  minimum  Richardson  number  in  the  flow  appears  to  be 
an  indicator  of  the  conditions  necessary  for  the  transfer  of  energy  to 
an  interfacial  disturbance.  In  other  words,  the  minimum  Richardson  num- 
ber appears  to  be  a good  indicator  of  the  degree  of  instability. 

Next,  an  explanation  of  why  the  high  frequency  waves  are  trapped  in 
the  interfacial  layer  is  given.  In  order  for  a disturbance  to  grow, 
the  Richardson  number  must  be  smaller  than  some  critical  value.  In 
addition,  only  disturbances  with  frequencies  less  than  the  natural  fre- 
quency (Brunt-Vai sala  frequency)  can  exist  in  the  fluid  as  pointed  out 
in  Chapter  1.  Thus  if  the  interfacial  layer  becomes  unstable  and  internal 
waves  of  certain  frequencies  (less  than  natural  frequency)  result,  these 
waves  are  confined  to  the  region  where  the  Brunt-Vaisala  frequency  is 
equal  to  or  greater  than  these  frequencies.  A definition  sketch  is  shown 
in  Figure  CO.  It  was  also  pointed  out  in  Chapter  1 that  disturbances 
propagating  in  a stratified  fluid  will  have  a vertical  component  of  pro- 
pagation until  it  reaches  a level  where  the  Brunt-Vaisala  frequency  coin- 
cides with  its  frequency.  Thus  disturbances  with  frequencies  less  than  N . 
(see  Figure  59)  will  propagate  into  the  stable  region  and  cease  to  exist. 

For  this  reason  only  waves  with  frequencies  greater  than  n . and  less  than 

m 1 n 

N are  believed  to  exist  in  the  interfacial  layer, 
max 

The  critical  value  of  the  Richardson  number  as  determined  by  these 
experiments  is  larger  than  that  predicted  by  the  linear,  inviscid  theory. 
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This  is  probably  duo  to  nonlinear  effects  caused  by  the  finite  size  of 

the  waves,  from  figure  59  the  critical  value  of  Richardson  number  (R.  ) 

' icr' 

is  approximately  10. 

The  changes  in  temperature  in  the  interfacial  region  were  small,  as 
can  be  seen  in  the  temperature  profiles  presented  in  Figure  61.  This 
coupled  with  the  fact  that  the  conductivity  readings  were  found  to  be 
only  slightly  dependent  on  temperature  (see  Figure  62)  made  temperature 
correction  of  the  conductivity  unnecessary. 

Unfortunately  the  attempt  to  analyze  the  frequency  spectrum  failed. 

The  frequencies  of  the  disturbances  were  expected  to  be  near  the  natural 
frequencies  in  the  interfacial  layer.  From  the  Brunt-Vaisala  frequencies 
plots  in  Figures  44  through  49,  it  can  be  seen  that  the  maximum  natural 
frequency  is  0.54  Hz.  From  the  visicorder  recording  of  velocity  (see 
Figure  73  and  74  in  the  Appendix),  the  frequency  of  the  disturbance  was 
estimated  to  be  between  0.5  and  0.67  Hz,  which  is  very  close  to  the  range 
of  natural  frequencies. 

In  these  experiments,  the  wave  steepness  was  not  excessively  high. 

It  is  expected  that  if  the  wave  steepness  is  much  larger,  the  disturbances 
will  be  in  the  turbulent  subrange  and  significant  changes  in  the  density 
profile  will  occur.  The  periodic  exchange  of  energy  between  the  internal 
and  interfacial  waves  should  no  longer  occur  when  this  happened. 
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CHAPTER  VII 

aspic  is  oi  i mum  studies 


This  thesis  only  investigated  the  initiation  of  the  high  frequency 
interfacial  waves  generated  by  long  internal  waves.  It  is  highly  rec- 
ommended that  the  interfacial  wave  spectrum  be  measured  and  analyzed 
in  order  to  determine  the  frequencies  that  exist  and  the  energy  level  of 
the  disturbances.  To  achieve  this,  a cylinder-type  hot-film  probe  is 
needed  for  measuring  very  low  fluid  velocities.  Also,  a low  noise  level 
I'M  tape  recorder  is  required. 

In  this  set  of  experiments,  complete  velocity  profiles  were  absent 
because  the  wedge-type  hot-film  probe  can  only  measure  the  velocity  in 
one  direction.  The  cylinder-type  hot-film  probe  will  also  allow  com- 
plete velocity  profiles  to  be  made. 

There  were  only  a limited  number  of  wave  frequencies  and  heights 
examined  in  tins  work.  This  should  be  expanded  to  include  more  internal 
wave  frequencies  and  a (treater  range  of  steepness  ratios.  With  suf- 
ficient data,  a critical  Richardson  number  will  be  established. 

The  growth  and  decay  of  the  internal  waves  should  be  studied  by 
measuring  the  spatial  as  well  as  the  temporal  changes  in  the  interfacial 
waves.  This  can  be  accomplished  by  moving  the  traverse  along  the  wave 
tank  and  by  monitoring  the  wave  heights  at  more  stations  using  the  wave 
gages . 

In  these  studies,  tho  transition  to  turbulence  is  not  examined.  It 
is  recommended  that  the  steepness  ratio  be  increased  and  the  transition 
vestigated.  Hopefully,  a transition  parameter  can  be  established. 
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